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Abstract 

We present updated predictions for the total production cross section of top-quark pairs at 
the Tevatron and at the LHC, and, at the LHC, of heavy-quark pairs with mass in the range 
0.5 — 2 TeV. For tt production at the LHC we also present results at = 10 TeV, in view of the 
expected accelerator conditions during the forthcoming 2008 run. Our results are accurate at the 
level of next-to-leading order in as, and of next-to- leading threshold logarithms (NLO-I-NLL). 
We adopt the most recent parametrizations of parton distribution functions, and compute the 
corresponding uncertainties. We study the dependence of the results on the top mass, and we 
assess the impact of missing higher-order corrections by independent variations of factorisation 
and renormalisation scales. 



*0n leave of absence from INFN, Sezione di Genova, Genoa, Italy. 
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1 Top production 



One of the first tasks of the Large Hadron Colhder (LHC) experiments when, later in 2008, they 
will start taking data, will be to re-discover the Standard Model. This will serve the double purpose 
of making sure that the detectors are well understood and work properly, as well as of improving 
the precision of previous measurements. With a predicted cross section for top quark pairs about 
a hundred times larger than at the Fermilab Tevatron, and a much higher design luminosity, the 
LHC is poised to become a real "top factory" pQ. This will allow for better measurements of the 
mass and the cross section. The former is expected to be measured with an ultimate uncertainty 
below 1 GeV [U [2] (to be compared with the most recent determination from the Tevatron, mt = 
172.6 lb 0.8 ± 1.1 GeV [3]). The cross section is expected to be measured within a year with a 15% 
accuracy, and eventually with an accuracy probably limited only by the knowledge of the LHC 
luminosity, expected to reach a precision of a few per cent [4] . 

Experimental measurements of the total cross section at the Tevatron [5l E] are usually compared to 
predictions [8', "9] compiled a few years agcQ. These predictions made use of the next-to-leading order 
(NLO) calculations of pjj, and of the soft-gluon next-to-leading-log (NLL) threshold resummation 
results obtained in [12^ [T3] and [13] respectively. Some logarithmic contributions of order higher 
than NLL have been included in the results of ref. [SI [15] and subsequent papers. In this way, while 
a complete NNLO calculation is still unavailable (the first ingredients, two loop virtual corrections 
in the ultrarelativistic limit [16^ I17j . fully massive results for the two-loop contribution in the 
quark-quark channel [181 US] and for the one- loop-squared contributions [201 IZI] , real emission at 
one loop [22], have been recently obtained), some NNLO contributions of order a|, of soft origin, 
can be obtained through an expansion and truncation of the Sudakov exponent. More recently, 
while this paper was being completed, the final ingredients required for a complete resummation 
of soft-gluon next-to-next-to-leading logarithms (NNLL) have been calculated [23j . and the relative 
impact on the total cross sections was explored. 

A key feature of our most recent study in [8] was an extensive exploration of the theoretical un- 
certainties affecting the prediction. The effect of the independent variations of renormalisation and 
of factorisation scales, which is the customary way to assess the impact of unknown higher-order 
contributions, was explored in detail. Parton distribution function sets (PDFs) providing a mean to 
estimate the associated uncertainty [23] were also used. It was determined in [8] that a significant 
fraction of the overall ±10 — 13% uncertainty in pp production at the Tevatron was originating 
from the PDFs, though higher orders were also contributing a fair share. This result should now be 
revisited on a number of counts. First, new PDF sets with errors, CTEQ6.5 [25], MRST2006nnlo 
|26j and CTEQ 6.6 [27] have appeared in the past few years. It is legitimate to wonder if they might 
come with a reduced uncertainty. Second, a similarly careful job of estimating the theoretical un- 
certainties for the best available prediction should be made for the LHC too. Third, since the most 
recent Tevatron measurements point to a lower mass than the central value mt = 175 GeV used in 
[8], it is useful to produce numerical predictions for an updated value of the top mass. Note that an 
analysis of the tt cross section has recently been performed in [27] . This study, carried out at the 
fixed-order, NLO level, focuses on the correlations of the top cross section with other observables, 
analyzed function of the PDF sets. 

^The precision of the cross section experimental measurements has recently become sufficiently good that extrac- 
tions of the top mass by comparing the measured cross section with the calculated value have become possible, and 
have been performed [10] , 
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We shall present our results in the form 

a = a(central)+t:! , (1) 

where u (central) is our best prediction, and Afi^-i- and AcpDp^ quantify the uncertainties due to 
higher perturbative orders and PDF choices, as specified in what follows. 

In order to streamline the calculation of the overall uncertainty (unknown higher orders and PDFs) 
we modify slightly the method employed in ref. [Sj and proceed as follows. 



Our best prediction u (central) is computed by setting the renormalisation and factorisation 
scales equal to m^, and with the central PDF set (within a given PDF error family). The 
cross section is calculated to NLO+NLL accuracy, exactly as in ref. [8]. 

The uncertainty on higher orders is estimated by varying the factorisation and the renormal- 
isation scales //p and //a independently around a central scale set by the top mass mt- We 
define the ratios 

^p = /ip/mt , = /iR/mj , (2) 

and we allow them to vary in the regions 0.5 < ^p,^r < 2, with the condition that 0.5 < 
?f/?r < 2. This means that none of the ratios ^ip/rrit, n^/nT't and /ip/^R can be larger than 
two or smaller than one-half, in order not to have in the perturbative expansion logarithms 
of arguments larger than a chosen (admittedly arbitrary) amount. Within this region the 
NLO-I-NLL cross section is evaluateql, and used to comput^ 



Ao"u+ = max 

{Cf.Cr} 



aiCF,^n)-cTihl)\, (3) 
a{^^,U)-a{l,l)]. (4) 



Ao"„_ = — min 

{Cf,5r} L 

All cross sections in these formulae are evaluated with the central PDF set (thus, a (1,1) = 
(T (central) here). We also introduce the symbols 

Scales-|- = o" (central) + A(T^+ (5) 
Scales— = (T(central) — A(t^_ (6) 

which we shall use in the following. 

By doing so we have established a variation interval of the cross section that can be considered 
as a reasonable estimate of the uncertainty due to unknown higher orders. It should be noted, 
however, that such an uncertainty should by no means be considered as distributed according 
to some probability law (for instance, with a Gaussian distribution) around the central value. 
In fact, it is more similar to a systematic than to a statistical uncertainty. This means that 
further arbitrary choices will have to be made in order to assign a 'confidence level' to this 
interval. 

Modern PDF sets come with a procedure to evaluate the propagation of their uncertainty 
onto a given physical observable. This is done by exploring the effect of using, along with a 



■^A NLL resummation function with independent renormalisation and factorisation scales is given explicitly in 
^The quantities A(t^+ and Aa"^_ are positive for all choices of top mass and scales we have considered. 
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'central' PDF set, a number of other sets (usually 40 for the CTEQ family PDFs, 30 for the 
MRST family ones) and properly combining their differences. According to the CTEQ and 
MRST Collaborations, the resulting uncertainty should roughly represent a 90% confidence 
level. We have chosen to follow the prescription by Nadolsky and Sullivan [29], and determine 
asymmetric uncertainties in the form 



Ac7pDF+ = J > I max 




a{set+i) - a{seto),a{set-i) - a{seto),0 ) , (7) 



AcTpoF- = J > I max 



-|\2 

a{setQ) — a{set^i),a{seto) — a{set^i),0 j . (8) 



where all cross sections are evaluated with 

^F = 1 , U = l. (9) 

In eqs. ([7]) and ([8]), seto represents the central set, and the sums run over all pairs of PDFs 
in the given PDF error set. For each pair, we denote by set+j and set^i the positive and 
negative displacement member of the pair. We also introduce the symbols 

PDFs+ = cr(central) + Ac7pdf+ (10) 
PDFs- = cr(central) - Acjpdf- (11) 

which we shall use in the following. 

To facilitate the determination of theoretical cross section corresponding to mass values different 
than the current best fit, we provide our results in the form of the coefficients of the parametrization 

cr{mt) =A + B{mt - 171) + C{mt - lllf + D{mt - lllf . (12) 

The parameters were fitted to the exact results in the range 150 <mt< 190 GeV, with a precision 
of the order of 1-2 per mille. The A coefficient has been fixed equal to the cross section at mt = 
171 GeV. The fit parameters are given in table [1] and table O for the Tevatron and the LHC 
respectivel}!!!. 

For each PDF set we have listed separately the 'central' value (scales = mt, central PDF set) and 
the maximum and the minimum found by varying the scales according to the above procedure and 
evaluating the asymmetric PDFs uncertainties. The main effects of the different PDFs and of the 
uncertainties can of course be read off directly from the A coefficient, which corresponds to the ti 
cross section evaluated at mt = 171 GeV. The display of results obtained with many different PDF 
sets, both very recent and older, is meant to allow for an easy estimate of the variation (or lack 
thereof) of the cross section predictions as a consequence of evolving parton distribution functions 
sets. 

We summarise here what might be considered our "best" predictions for ti production at the LHC, 
at mt = 171 GeV: 

^NLO+NLL(LHC,m4 = 171 GeV,CTEQ6.5) = 908 tf^%f%) (scales) ^sglSo) (PDFs) pb (13) 



■'Since the LHC is scheduled to run in 2008 at a centre of mass energy of 10 TeV, we have also provided predictions 
for this energy in table O in the case a measurement of the total cross section of top production should prove possible. 
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Xevatron, pp at a 


/i = I960 GeV 


A fnb) 


B (pb/GeV) 


C (pb/GeV^) 


D (pb/GeV^) 




Central 


7.59 


-0.237 


4.39 xlO"^ 


-6.32 xlO"^ 




Scales+ 


7.89 


-0.247 


4.60 xlO"^ 


-6.66 xlO'^ 


CTEQ6M 


Scales — 


7.07 


-0.221 


4.11 xlO"^ 


-5.92 xlO~^ 




PDFs+ 


8.26 


-0.260 


4.86 xlO"^ 


-7.02 xlO"^ 




PDFs- 


7.12 


-0.222 


4.08 xlO"^ 


-5.82 xlO"^ 




Central 


7.77 


-0.244 


4.53 xlO^^ 


-6.51 xlO"^ 




Scales+ 


8.08 


-0.254 


4.74 xlO"^ 


-6.86 xlO"^ 


CTEQ6.1 


Scales — 


7.23 


-0.227 


4.23 xlO"^ 


-6.09 xlO"^ 




PDFs+ 


8.53 


-0.269 


5.04 xlO"^ 


-7.27 xlO"^ 




PDFs- 


7.20 


-0.224 


4.12 xlO"^ 


-5.87 xlO~^ 




Central 


7.61 


-0.237 


4.38 xlO"^ 


-6.28 xlO"^ 




Scales+ 


7.90 


-0.247 


4.58 xlO"^ 


-6.61 xlO"^ 


CTEQ6.5 


Scales — 


7.08 


-0.221 


4.10 xlO"^ 


-5.89 xlO'^ 




PDFs+ 


8.14 


-0.256 


4.78 xlO"^ 


-6.91 xlO"^ 




PDFs- 


7.24 


-0.224 


4.11 xlO"^ 


-5.85 xlO~^ 




Central 


7.48 


-0.233 


4.32 xlO"^ 


-6.20 xlO"^ 




Scales + 


7.77 


-0.243 


4.52 xlO-2 


-6.53 xlO"^ 


CTEQ6.6 


Scales — 


6.96 


-0.218 


4.04 xlO"^ 


-5.80 xlO~^ 




PDFs + 


7.99 


-0.251 


4.70 xl0"3 


-6.79 xlO~^ 




PDFs - 


7.09 


-0.220 


4.02 xlO"^ 


-5.72 xlO"^ 




Central 


7.66 


-0.242 


4.53 xlO"^ 


-6.60 xlO~^ 




Scales+ 


7.97 


-0.252 


4.75 xl0~2 


-6.98 xlO"^ 


MRST2001E 


Scales— 


7.13 


-0.225 


4.24 xlO-3 


-6.17 xlO"^ 




PDFs+ 


7.94 


-0.252 


4.75 xlO"^ 


-6.95 xlO"^ 




PDFs- 


7.44 


-0.233 


4.35 xlO"^ 


-6.31 xlO"^ 


MRST2004nlo 


Central 


7.99 


-0.253 


4.77 xlO"^ 


-6.95 xlO"^ 




Central 


7.93 


-0.253 


4.76 xlO"^ 


-6.92 xlO-^ 




Scales+ 


8.27 


-0.264 


5.00 xlQ-^ 


-7.33 xlO-5 


MRST2006nnlo 


Scales— 


7.37 


-0.235 


4.44 xlO"^ 


-6.45 xlO"^ 




PDFs+ 


8.17 


-0.261 


4.93 xl0^3 


-7.19 xl0~5 




PDFs- 


7.73 


-0.245 


4.61 xlO-3 


-6.68 xlO"^ 



Table 1: Coefficients of the parametrization ^(mf) = A-\-B{mt — ni) + C{mt — 171)^ + D{mt — 171)^ 
for tlie NLO-I-NLL tt cross section (picobarn) at the Tevatron, for various PDF sets. The fit must 
not be used outside the range 150 < mt < 190 GeV. The quantities Scalesit and PDFsit are defined 
in eqs. (0), ©, ([ID]), and (HH). 
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LHC, pp at ^J~s - 


= 10 TeV 


A (pb) 


B (pb/GeV) 


C (pb/GeV^) 


D (pb/GeV^) 




Central 


425 


-12.1 


0.211 


-2.89 xlO^^ 




Scales+ 


462 


-13.2 


0.232 


-3.20 xlO"^ 


CTEQ6M 


Scales — 


386 


-10.9 


0.189 


-2.58 xlO"^ 




PDFs+ 


445 


-12.5 


0.216 


-2.94 xlO"^ 




PDFs- 


406 


-11.7 


0.205 


-2.82 xlO"^ 




Central 


428 


-12.1 


0.211 


-2.87 xlO"^ 




Scales -|- 


465 


-13.2 


0.232 


-3.19 xlO"''^ 


CTEQ6.1 


Scales — 


389 


-10.9 


0.189 


-2.57 xlO"^ 




PDFs+ 


450 


-12.5 


0.216 


-2.93 xlO"^ 




PDFs- 


406 


-11.6 


0.205 


-2.81 xlO~^ 




Central 


414 


-11.7 


0.205 


-2.79 xlO"^ 




Scales -|- 


450 


-12.9 


0.226 


-3.09 xlO"^ 


CTEQ6.5 


Scales — 


376 


-10.6 


0.184 


-2.50 xlO"^ 




PDFs+ 


434 


-12.2 


0.211 


-2.85 xlO"^ 




PDFs- 


396 


-11.3 


0.199 


-2.72 xlO"^ 




Central 


414 


-11.8 


0.206 


-2.81 xlO"^ 




Scales + 


451 


-12.9 


0.227 


-3.12 xlQ-^ 


CTEQ6.6 


Scales — 


376 


-10.6 


0.185 


-2.51 xlO"^ 




PDFs + 


433 


-12.2 


0.211 


-2.86 xlO"^ 




PDFs - 


396 


-11.4 


0.200 


-2.75 xlO"''^ 




Central 


446 


-12.6 


0.217 


-2.94 xlO"^ 




Scales+ 


486 


-13.8 


0.240 


-3.27 xlQ-^ 


MRST2001E 


Scales— 


405 


-11.3 


0.195 


-2.63 xlQ-^ 




PDFs+ 


457 


-12.8 


0.220 


-2.97 xlQ-"' 




PDFs- 


439 


-12.4 


0.216 


-2.92 xlQ-^ 


MRST2004nlo 


Central 


455 


-12.8 


0.221 


-2.99 xlO"^ 




Central 


446 


-12.5 


0.216 


-2.92 xlO"^ 




Scales+ 


486 


-13.7 


0.238 


-3.24 xlQ-^ 


MRST2006nnlo 


Scales— 


404 


-11.3 


0.194 


-2.60 xlQ-^ 




PDFs+ 


454 


-12.7 


0.218 


-2.93 xlO-3 




PDFs- 


438 


-12.3 


0.214 


-2.89 xlQ-^ 



Table 2: Coefficients of the parametrization a{mt) = A+B{mt-ni) + Cimt-nif + D (mf — 171) 
for the NLO+NLL tt cross section (picobarn) at the LHC with ^/s = 10 TeV, for various PDF 
sets. The fit must not be used outside the range 150 < mt < 190 GeV. The quantities Scalesit and 
PDFsib are defined in eqs. dSj), (fTO]) . and ([TT]) . 
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LHC, pp at ^J~s - 


= 14 TeV 


A (pb) 


B (pb/GeV) 


C (pb/GeV^) 


D (pb/GeV^) 




Central 


933 


-25.3 


0.423 


-5.60 xlO^^ 




Scales+ 


1018 


-27.7 


0.468 


-6.22 xlO"^ 


CTEQ6M 


Scales — 


846 


-22.8 


0.379 


-4.99 xlO"^ 




PDFs+ 


962 


-25.8 


0.432 


-5.73 xlO"^ 




PDFs- 


903 


-24.6 


0.413 


-5.44 xlO"^ 




Central 


934 


-25.2 


0.421 


-5.56 xlO"^ 




Scales+ 


1019 


-27.7 


0.466 


-6.19 xlO"''^ 


CTEQ6.1 


Scales — 


847 


-22.7 


0.377 


-4.95 xlO"^ 




PDFs+ 


965 


-25.8 


0.430 


-5.70 xlO"^ 




PDFs- 


902 


-24.5 


0.411 


-5.40 xlO"^ 




Central 


908 


-24.5 


0.411 


-5.46 xlO"^ 




Scales -|- 


990 


-26.9 


0.455 


-6.08 xlO"^ 


CTEQ6.5 


Scales — 


823 


-22.1 


0.368 


-4.87 xlO"^ 




PDFs+ 


938 


-25.2 


0.420 


-5.57 xlO"^ 




PDFs- 


879 


-23.9 


0.401 


-5.29 xlO"^ 




Central 


911 


-24.7 


0.413 


-5.47 xlO~^ 




Scales + 


993 


-27.1 


0.457 


-6.09 xlQ-^ 


CTEQ6.6 


Scales — 


826 


-22.2 


0.370 


-4.87 xlO"^ 




PDFs + 


939 


-25.2 


0.422 


-5.58 xlO"^ 




PDFs - 


881 


-24.0 


0.404 


-5.36 xlO"''^ 




Central 


965 


-25.9 


0.429 


-5.63 xlO"^ 




Scales+ 


1054 


-28.4 


0.475 


-6.27 xlQ-^ 


MRST2001E 


Scales— 


874 


-23.3 


0.384 


-5.00 xlQ-^ 




PDFs+ 


981 


-26.2 


0.434 


-5.68 xlQ-"' 




PDFs- 


954 


-25.6 


0.426 


-5.57 xlQ-^ 


MRST2004nlo 


Central 


982 


-26.3 


0.436 


-5.72 xlO"^ 




Central 


961 


-25.7 


0.426 


-5.58 xlO"^ 




Scales+ 


1050 


-28.3 


0.472 


-6.21 xlQ-^ 


MRST2006nnlo 


Scales— 


870 


-23.1 


0.381 


-4.96 xlQ-^ 




PDFs+ 


972 


-25.9 


0.428 


-5.62 xlQ-^ 




PDFs- 


949 


-25.4 


0.422 


-5.53 xlQ-^ 



Table 3: Coefficients of the parametrization (T(mt) = A-\-B{mt — 171) + C{mt — 171)^ + D{mt — 171)^ 
for the NLO+NLL tt cross section (picobarn) at the LHC, for various PDF sets. The fit must not 
be used outside the range 150 < mt < 190 GeV. The quantities Scalesit and PDFsit are defined in 
eqs. ©, (©, (Uni), and (HH). 
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^NLO+NLL^LHC,mf = 171 GeV, MRST2006nnlo) = 961 ^gilgj^j (scales) ^JsSlS) (PDFs) pb 

(14) 

Note that we quote separately the results for MRST and CTEQ, since they are not fully consis- 
tent. For reference, we also give here the pure fixed-order results (i.e., without including threshold 
resummation) at the NLO and the LO 

a,^-LO(LHC,mt = 171 GeV,CTEQ6.5) = 875 ^JJoSIlS) (^^^^^s) -29(3:3%) (P^^^) (1^) 

(Tjf°(LHC,mt = 171 GeV,CTEQ6.5) = 583 tlfo^l'jyj) (scales) t^^^^iyj) (PDFs) pb (16) 



af^^^{LRC,mt = 171 GeV, MRST2006nnlo) = 927 ^}J7[}}:5%J (scales) ^}2(L3%) (PDFs) pb (17) 

a^p{LRC,mt = 171 GeV, MRST2006nnlo) = 616 tS^Hy^j (scales) ^7 8(} 3%) (PDFs) pb (18) 

We have decided not to combine the scales and PDFs uncertainties into a single error. The reason 
for refraining from doing so is that the scales uncertainty (and, to some extent, probably also the 
PDFs one) is not fully characterised in statistical terms. As a consequence, additional hypotheses 
will be needed in order to combine the two uncertainties into a single probability density function for 
the resulting cross section, with well defined confidence levels. Further discussions on the interplay 
between scales and PDFs uncertainties can be found in Appendix lAl 

We finally present our "best" predictions for tt production at the Tevatron, at nit = 171 GeV: 

°+''^'^(Tev,mt = 171 GeV,CTEQ6.5) = 7.61 tH^'^J) (scales) (PDFs) pb (19) 

^NLO+NLL(rpg^^^^ = 171 GeV,MRST2006nnlo) = 7.93 ^0 56(7?5j (scales) ^.olf^^^J) (PDFs) pb . 

(20) 



As done for the LHC in eqs. (|T5|) - (|T8|) . we also report the NLO and LO results: 

4LO(Tev,mt = 171 GeV,CTEQ6.5) = 7.35 (scales) tHZSJ) (PDFs) pb (21) 

4°(Tev,mt = 171 GeV,CTEQ6.5) = 5.92 !?:Jlg^:;5J (scales) tlfJ^^'S^ (PDFs) pb (22) 



aj^-^"(Tev,mt = 171 GeV, MRST2006nnlo) = 7.62 _o:88(n 6%) (scales) _o:i8(2.4%) (PDFs) pb (23) 
a'^°{Tev,mt = 171 GeV, MRST2006nnlo) = 6.05 ti^i^S^^^J^ (scales) ^o;}3|2:55) (PDFs) pb . (24) 

2 Discussion of the ti cross section results at the LHC 

For all parameter choices we have considered, the scales uncertainties affecting the tt cross section at 
the LHC are much larger than those due to PDFs. In this section, we therefore focus on exploring 
the effect of the NLL resummation on the cross section. We can do so by comparing the NLO 
prediction at = 171 GeV, and its uncertainty due to scales variations as described above, in the 
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NLO and NLO+NLL approximations respectively. We find that the 'central' prediction is increased 
by less than 4%, and the scales uncertainty is only very mildly affected, going from ±11.5% in the 
NLO case to ±9% in the NLO+NLL one. This points to a relatively minor impact of threshold 
resummation on the LHC cross section, as expected as a consequence of the relatively large distance 
of the tt production threshold from the LHC centre of mass energy (for comparison, the uncertainty 
at the Tevatron is almost halved when going from NLO to NLO+NLL). One should also note that, 
again contrary to the Tevatron case, exploring independent scale variations has a non-negligible 
effect: keeping /Ur = /Xp (as done e.g. in ref. [27]) would result in an uncertainty estimate for the 
NLO+NLL case of only We remind the reader that the fact PDF fits are performed with 

/iR = /Xp does not force us to use //r = //p in the cross section. In fact an independent variation 
of the two scales in our matched calculation leads to variations in the result that are beyond the 
NLO+NLL approximation. It is thus legitimate to add this independent variation to the sources 
of uncertainties. It then turns out that the //r ^ //p approach leads to a much larger variation. 
We thus conclude that there may be accidental cancellation in the scale variation when one keeps 
/iR = /ip, leading to an unreliably small estimate of the error. 

Another important element in the assessment of the systematics related to the resummation is the 
estimate of the impact of beyond-NLL corrections. To parametrize these corrections, a constant A 
was introduced in [13J (where more details about its role are given): 

Cij^Cij (i- . ij = qq,99- (25) 

Cij here represents the A^-independent term of the Mellin transform of the NLO partonic cross 
section. The replacement in the previous equation gives vanishing first moments, is irrelevant for 
large A^, and does not introduce poles on the real N axis. Different choices of A give rise to different 
resummed cross sections, all consistent with each other at the NLL level and NNLL leveH. They 
therefore parametrize the possible exponentiation of finite, non-logarithmic terms appearing at 
orders higher than NLO. It was noticed already in [13] that the choice ^ = was leading to a scale 
dependence typically a factor of two smaller than was obtained with A = 2. For the sake of being 
conservative, we therefore selected A = 2 in our subsequent phenomenological analysis [8j, as well 
as in the results presented in the previous section. We would like to reiterate here this observation, 
by showing the results that we would have obtained at the LHC if we had chosen A = 0. In the 
case of the CTEQ6.5 PDF, and mt = 171 GeV, the scale dependence obtained by varying fip and 
/lxr independently is: 

^NLO+NLL(A=o)^LHC, mi = 171 GeV, CTEQ6.5) = 945 ^85(9°o%) (scales) pb [/ip ^ //rJ , (26) 

which is approximately 5% larger than the value obtained with A = 2, a variation consistent with 
the estimated uncertainty. If we had chosen to set fip = /Xr, the result would have been: 

^NLO+NLL(A=o)^-^jj(^^^^ = 171 GeV, CTEQ6.5) = 945 ^J^So) (=^^^1^^) [-"^ = ■ (^7) 

Notice that the combination of A = and /Xr = /Xp leads to a dramatic reduction of the uncertainty. 
In particular, the reduction is significant also with respect to the A = 2, /Xr = /Xp case, discussed 
above. The choice ^4 = and /Xr = /Xp is what was used in the recent analysis of the resummed 



^Changing A corresponds to vary terms suppressed by powers of A*', which cannot be determined within the soft 
;luon approximation. 
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NNLL cross section of |23j , leading to a similar uncertainty, at the 2% level. We conclude that it is 
not as yet clear whether the improvement found in [23j is a genuine reduction of the uncertainty, 
that would survive the independent variation of and /ip and with the introduction of a parameter 
similar to our A. 

We also note that the resummed NLO+NLL results quoted in ref. [23] differ from ours to the extent 
of a few percent. We have checked that this is not due to the choice made in [23J of limiting the 
contribution of NLL resummation to a tt mass range near the production threshold, and using only 
the NLO result above such range (contrary to what stated in [23], we do not perform the matching 
in this way, but follow instead the procedure detailed in [l3]). Rather, the small discrepancy is 
due to the misleading way in which eq. (125p is presented in [13], and in which it was accordingly 
interpreted and used in ^23]. There is in fact a mismatch in the Mellin N argument appearing in 
the expression of coefficient Cij used in eqs. (54) and (58) of [13], with + 1 being the correct 
argument of the resummation function, rather than N. The numerical implementation of these 
relations, in [13] and in this paper, are nevertheless consistent, and equivalent to a shift ^ + 1 
in (I25p . An erratum to clarify this issue has been submitted for ^13j. 

Another issue we wish to comment on is the relative size of the scale and the PDFs uncertainty. 
It was observed in [8] that the latter was important, and almost dominant, at the Tevatron. This 
appears not to be the case anymore at the LHC: again according to the procedure described above, 
at nit = 171 GeV we find uncertainties of the order of ±3% with CTEQ6.5 and ±1.5% with 
MRST2006nnlo and MRST2001E. The main reason for this improvement is that at the LHC the 
range of x values for the partons relevant to top production is much smaller than at the Tevatron, 
and falls in a region where the experimental knowledge of both quark and gluon PDFs is much 
better constrained by data. It is worth noting that the central results given by these two PDF sets, 
differing by about 6%, are not fully compatible, despite (or because of) the apparently very small 
estimated uncertainty. We also point out the MRST and CTEQ use different conventions for the 
Tolerance parameter; the consequence of this is that, had the two collaborations followed exactly 
the same fitting procedure, the PDF uncertainty resulting from using an MRST family set would 
still be a factor of about \/2 smaller than that obtained with a CTEQ set. 

While the PDFs uncertainty is probably still somewhat underestimated, as shown by the partially 
conflicting central values of CTEQ6.5 and MRST2006nnlo, it is probably safe to conclude that, the 
very interesting progress with the NNLL resummation [23j notwithstanding, a deflnitive assessment 
of our understanding of the tt cross section at the LHC will have to wait for the full, massive NNLO 
calculation. 

3 Very heavy fermion production 

We now present production rates for a pair of fermions (belonging to the fundamental representation 
of SU {2>) colour) heavier than top, using the same computations described for the cross sections. 
We shall generically denote such fermion pair by TT. These particles arise naturally in BSM 
theories with strongly-coupled dynamics; they can be of different species, which can for example be 
classified according to their transformation properties under SU(2)Li8'SU(2)R(g)U(l). As far as pair 
production is concerned, however, these details are largely irrelevant, since this process is expected 
to be dominated by QCD effects, and this is the reason why we can apply our NLO+NLL, NLO, 
and LO results to the computations of TT rates. In doing so, we shall neglect possible contributions 
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rrix 


NLO+NLL 


NLO 


LO 


0.5 


+232(5.8%) +466(11.7%) 
40U6_27g(g -332(8.3%) 


+342(9%) +455(12.0%) 
'^oU2_421(ll.l%) -322(8.5%) 


„„r,£;+876(32.1%o) +314(11. 5%o) 
^'^"-618(22.7%) -221(8.1%) 


0.6 


.r,Q+76.3(5.3%) +195(13.7%) 
^^^^-93.2(6.5%) -134(9.4%) 


, „^^+116(8 6%il +188(14 0%) 
l'30^_148(ii.o%) -129(9.5%) 


^+319(32 5%) +130(13 3%) 
ySU./_225(22.9%) -88.7(9.1%) 


0.7 


c-yy ^+28.7(5%) +89.0(15.4%) 
^''•"-36.0(6.2%) -59.6(10.3%) 


+44.9(8.2%) +85.5(15.7%)) 
^^^•^-59.3(10.9%) -57.0(10.5% 


+131(32.9%) +59.0(14.8%) 
•399.i_g2 ^23.1%) -39.2 9.8%) 


0.8 


^+12.0(4.7%) +43.4(17%) 
^"^"•^-15.4(6%) -28.4(11.1%) 


^+19 3 8%o) +414 17 2%)) 
^*^-^-26.1(10.8%) -27.0(11.2%) 


„+58. 9(33.1%) +28.5(16.1%) 
1' '•0-41.3(23.2%) -18.6(10.5%) 


0.9 


7+5.41(4.4^ 
^^-^•'-7.06(5.8^ 


7c) +22.3(18.4%) 
ii) -14.4(11.8%) 


+8.95(7.87„j +21.2(18.5%)) 
114.3_-^2.3(10.8%) -13.5(11.9% 


+28.3(33.4%)) +14.6(17.2%)) 
o4.71_-^gg^23.4%) -9.32(11%) 


1.0 


,„+2. 59(4.2%) +12.0(19.6%) 
"^•-■^^-3.45(5.6%) -7.58(12.4%) 


^„„^+4 42 7 7%) +113 19 7%)) 
^'•^^-6.17(10.8%) -7.11(12.4%) 


^„+14.3(33.7%) +7.75(18.2%) 
^^•^'-10.0(23.5%) -4.88(11.5%) 


1.1 


or, pci+l. 32(4.1%) +6.68(20.8%) 
"^^•^^-1.76(5.5%) -4.15(13%) 


.+2.29(7.7%)) +6.25(20.9%) 
29. y4__3 24(10.8%) -3.88(12.9%) 


+7.56(33.9%,) +4.28(19.2%)) 
^^•■^^-5.28(23.6%) -2.66(11.9%) 


1.2 


„ . +0.706(4.1%) +3.83(22%) 
^'•'*^-0.939(5.4%) -2.35(13.5%) 


„ „„+l 24(7 6%) +3 57(22%) 
^"•^'^-1.76(10.9%) -2.18(13.4%) 


+4.13(34.2%)) +2. 43(20. 17c) 
12.10_2.88(23.8%) -1.50(12.4%) 


1.3 


„ „„„+0.388(4%) +2.25(23.2%) 
^•'"^'-0.516(5.3%) -1.36(14%) 


+0.693(7.7%)) +2.09(23.1%)) 
^■^^^-0.989(10.9%) -1.26(13.9%) 


-+2.32(34.4%) +1.42(21.1%) 
"•'^^-1.61(23.9%) -0.864(12.8%) 


1.4 


c.„„+0.218(3.9%) +1.36(24.3%) 
D.D ' o_o.291(5.2%) -0.810(14.5%) 


+0.398(7.7%) +1.25(24.2%) 
^•^"^-0.569(11%) -0.745(14.4%) 


+1.34(34.7%) +0.850(22.1%) 
3.848_g g27(24.i%) -0.511(13.3%) 


1.5 


o r,f,f.+0.126(3.9%) +0.833(25.5%) 
•^•^"^-0.168(5.2%) -0.492(15.1%) 


„ „+0 235(7 8%) +0 763(25 3%) 
"^•^^-^-0.335(11.1%) -0.450(14.9%) 


„ „„„+0 783(35%)) +0 518(23 1%)) 
^•■^■^O-0. 543(24.2%) -0.309(13.8%) 


1.6 


1 000+0.074(3.8%) +0.520(26.8%) 
^•^"^°-0.099(5.1%) -0.304(15.7%) 


„op+0. 141(7.9%) +0.474(26.5%) 
1.785_Q 200(11.2%) -0.277(15.5%) 


+0.467(35.3%) +0.321(24.3%,) 
1.323_g 323(24.4%) -0.190(14.4%) 


1.7 


1 1 «n+0-044(3.7%) +0.329(28.2%,) 
^•-■^"^-0.059(5.1%) -0.191(16.3%) 


„ +0.086(8%) +0.299(27.8%)) 
i.U ' '3_o.l22(11.4%) -0.173(16.1%) 


_ +0.282(35.6%)) +0.202(25.5%)) 
^•'^■^-0. 195(24.6%) -0.119(15%) 


1.8 


.+0.026(3.7%) +0.212(29.6%) 
^•'^^-0.036(5%) -0.123(17.2%) 


„ „j.„+0.053(8.2%) +0.191(29.2%) 
^•"^"^-0.075(11.5%) -0.109(16.8%) 


n /fon+0-173(35.9%) +0.129(26.9%) 
^•'^o^-0. 119(24.7%) -0.075(15.7%) 


1.9 


440+0-016(3.7^ 
^•^^^-0.022(5% 


^} +0.137(31.2%) 
-0.078(17.7%) 


n zLni +0-033(8.4%) +0.123(30.8%) 
^•^^^-0.047(11. 7%) -0.070 17.5%) 


n 904+0-107(36.3%) +0.083(28.3%) 
'-'•^^^-0.073(24.9%) -0.048(16.4%) 


2.0 


n 974+0.010(3.6^ 
^•^'^-0.013(5% 


f'o) +0.090(32.9%) 
-0.051(18.6%) 


n 940+0.021(8.5%) +0.080(32.4%) 
U.Z48_Q 029(11.8%) -0.045(18.3%) 


p. , 0-1+0.066(36.6%) +0.054(30%) 
^•-'-^-'--0. 045(25.1%) -0.031(17.3%) 



Table 4: Cross sections (in fb) for the production of TT pairs at the LHC, computed with CTEQ6.5 
PDFs. The mass of the heavy fermion T is expressed in TeV. For each entry of the table, we give 
the central value of the cross section, with the scale and PDF uncertainties. 



of non-SM intermediate states resulting from e.g. a qq annihilation. Our aim is therefore not that 
of providing a complete phenomenological study of TT cross sections at the LHC, but rather that 
of assessing the scale and PDF uncertainties affecting the QCD contribution to the production of 
heavy fermion pairs. In what follow, we shall consider the mass range 0.5 TeV< rriT < 2 TeV for 
the heavy fermion. 

Our results are presented in tables [Hand [5l A few comments are in order. 

• The scale dependence of the NLO+NLL cross section is small for all values of rriT, and 
decreases monotonically with increasing niT- This is to be expected, since the heavier the 
fermion, the closer the kinematics is to the threshold. 

• The scale dependence of the NLO cross section starts by decreasing, but then tends to increase 
with increasing mx- This is the signal of the necessity of including threshold corrections. On 
the other hand, the scale dependence of the LO cross section is always extremely large. This 
fact must be taken into proper account if an estimate of the K factor is needed. In particular, 
one must precisely understand which hard scale is used in a LO computation (e.g. in a 
standard parton-shower Monte Carlo). 

• The relative PDF uncertainty is extremely large in the case of CTEQ6.5. When MRST2006nnlo 
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rriT 


NLO+NLL 


NLO 


LO 


0.5 


,^^„+267.4(6%) +197.4(4.4%) 
^^"^-314.0(7%) -172.6(3.9%) 


.90f;+392.9(9.3%) +191.9(4.5%) 
^"^^^ -480.5^^-3%) -167.7(4%) 


or,, „+988.3(32.8%) +132.0(4.4%) 
"^^-^'-694.6(23%) -115.9(3.8%) 


0.6 


1 c:qq+88.7(5.5%) +81.2(5.1%) 
-'-Oy»_io7.0(6.7%) -70.1(4.4%) 


-. 0+134.9(8.9%) +78.2(5.2%) 
l-Ol-'3_170.9(ii.3%) -67.5(4.5%) 


r.oQ+362.9(33.3%) +53.8(4.9%) 
lUOy_253 g(23,3%) -46.6(4.3%) 


0.7 


8+33.7(5.2%) +36.3(5.6%) 
O4o-o_41.6(6.4%) -31.0(4.8%) 


, (j+52.9(8.6%) +,34.7(5.7%) 
01-1-»_68.9(11.3%) -29.6(4.8%) 


AAA 0+1^5.0(33.8%) +23.9(5.4%) 
^^^•^-104.5(23.5%) -20.4(4.6%) 


0.8 


9S« 9+14.1(4.9%) +17.3(6%) 
^oo-^_17.8(6.2%) -14.6(5.1%) 


P^+22.9(8.4%) +16.4(6.1%) 
^'1-U_30.5(11.2%) -13.9(5.1%) 


f^+67.5(34.1%) +11.4(5.7%) 
iyo-l^_46.9(23.7%) -9.57(4.8%) 


0.9 


, „„ „+6.42(4.7%) +8.72(6.4%) 
-■-^'•^-8.21(6%) -7.28(5.3%) 


1 9« £;+10.7(8.3%) +8.23(6.4%) 
I^o-O_14.5(ll,3%) -6.85(5.3%) 


QA 04+32.5(34.4%) +5.74(6.1%) 
y4.J4_22 5(23.9%) -4.72(5%) 


1.0 


fiS 07+3.09(4.5%) +4.61(6.7%) 
OO-y'_4.02(5.8%) -3.77(5.5%) 


fiA /1 0+5.30(8.2%) +4.32(6.7%) 
D4.4»_.^27(11.3%) -3.52(5.5%) 


.0+16.5(34.8%) +3.05(6.4%) 
^'■^^-11.4(24%) -2.44(5.1%) 


1.1 


0^.91+1.56(4.3%) +2.53(7%) 
^D./l_2.o6(5.7%) -2.03(5.6%) 


00 7r+2.76(8.2%) +2.37(7%,) 
05.^0_3_g2(ii,3%) -1.88(5.6%) 


9zL 07+8-72(35%) +1.69(6.8%) 
^^■O' -6.02(24.2%) -1.31(5.3%) 


1.2 


1 Q fiQ+0-824(4.2%) +1.44(7.3%) 
-^^•"^-1.10(5.6%) -1.12(5.7%) 


,g on+l-49(8.2%) +1.35(7.4%) 
-^°-^^-2.08(11.4%) -1.04(5.7%) 


1 q rf.+4.77(35.3%) +0.987(7.3%) 
i>i.i3U_3 29(24.4%) -0.734(5.4%) 


1.3 


^0+0.449(4.1%) +0.857(7.8%) 
-'--l-Uo_Q go4(5 5%) -0.642(5.8%) 


99+0.837(8.2%) +0.802(7.9%) 
-^^•^^-1.17(11.5%) -0.593(5.8%) 


„ (-00+2.68(35.6%) +0.597(7.9%) 
'•0'3^-1.85(24.5%) -0.425(5.6%) 


1.4 


a oo/i +0.252(4%) +0.524(8.3%) 
"•^^^-0.341(5.4%) -0.377(6%) 


c 0+0.482(8.3%) +0.492(8.4%) 
0-o4'3_0.674(11.5%) -0.348(6%) 


. onq+l-55(35.9%) +0.372(8.7%) 
4-^U^-1.06(24.7%) -0.253(5.9%) 


1.5 


7^,9+0.145(3.9%) +0.329(8.9%) 
^•'^^-0.197(5.3%) -0.227(6.1%) 


.p>Q+0.284(8.3%) +0.310(9.1%) 
"^■^^^-0.397(11.6%) -0.210(6.2%) 


9 rr|f;+0-908(36.2%) +0.238(9.5%) 
^.OUO_gg22(24.8%) -0.155(6.2%) 


1.6 


9 9nQ+0-085(3.8%) +0.211(9.6%) 
^•^U^_0.ll6(5.3%) -0.140(6.3%) 


9 n99+0-l'''l(8.5%) +0.199(9.9%) 
Z.U/Z_Q 238(11.8%) -0.130(6.4%) 


1 /ioq+0.542(36.5%) +0.155(10.4%) 
l-4o>J-0.370(25%) -0.097(6.5%) 


1.7 


1 qon+0-051(3.8%) +0.138(10.4%) 
-'-■^''^-0.069(5.2%) -0.088(6.6%) 


. 9,„+0.104(8.6%) +0.131(10.7%) 
1-^1 ' -0.145(11.9%) -0.082(6.7%) 


n oQn+0-328(36.8%) +0.102(11.5%) 
U.»yU_Q_224(25.1%) -0.061(6.9%) 


1.8 


p, 0+0.031(3.8%) +0.091(11.2%) 
U.»i^_Q Q42(5.2%) -0.056(6.9%) 


„ +0.065(8.8%) +0.087(11.7%) 
^•'^-'--0.089(12.1%) -0.052(7%) 


n rACt+^-'^^^i^'^-'^"'^'') +0.068(12.6%) 
^•^^^-0.137(25.3%) -0.039(7.3%) 


1.9 


n c:n9+0-019(3.7%) +0.061(12.2%) 
U.OU/_Q Q26(5.1%) -0.036(7.2%) 


„ .(.(.+0.041(8.9%) +0.058(12.8%) 
^•^^^^-0.056(12.2%) -0.034(7.4%) 


n qqn+0-124(37.5%) +0.046(13.8%) 
^•"^^^-0.084(25.5%) -0.026(7.8%) 


2.0 


r, 019+0-012(3.7%) +0.041(13.2%) 
^•"^-^^-0.016(5.1%) -0.023(7.5%) 


n 909+0-026(9.1%) +0.039(13.9%) 
^•^°^-0.035(12.4%) -0.022(7.8%) 


n 9n/l+0-077(37.9%) +0.031(15%) 
^•^^^-0.052(25.7%) -0.017(8.2%) 



Table 5: Cross sections (in fb) for the production of TT pairs at the LHC, computed -with 
MRST2006nnlo PDFs. The mass of the heavy fermion T is expressed in TeV. For each entry 
of the table, we give the central value of the cross section, with the scale and PDF uncertainties. 
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sets are used, that uncertainty is smaller by a factor of about 2-3, consistently with what 
already observed in the case of top production. By and large, the PDF uncertainty affects 
equally the NLO+NLL, the NLO, and the LO cross sections. At the largest my values consid- 
ered here, it prevents one from giving a precise prediction even in the case of the NLO+NLL 
computation. It will therefore be highly desirable to measure the PDFs at the LHC for 
intermediate- and large-x regions, using e.g. low-mass final states produced at large rapidity. 

4 Conclusions 

In this paper we have produced and tabulated updated predictions for the next-to-leading order 
plus next-to-leading log resummed (NLO+NLL) cross sections for ti production at the Tevatron 
and at the LHC. QCD cross sections for heavy fermion production at the LHC are also given. 

The theoretical uncertainties due to unknown higher orders and to the imperfect knowledge of the 
parton distribution function sets are explored in detail and also tabulated. NLO and LO results 
are also given, for reference and comparison. 

The main results for ti production at the LHC (v^ = 14 TeV) are the following: 

^NLO+NLL (LHC, = 171 GeV,CTEQ6.5) = 908 tlf^'S!) (scales) :!:29(3:2%) (PDFs) pb (28) 

^NLO+NLL (LHC, mt = 171 GeV, MRST2006nnlo) = 961 ^gijgj^j (scales) ^JaSl'S) (PDFs) pb 

(29) 

Cross sections obtained with two of the most recent PDF sets are given, since they appear to be 
only partially compatible within their respective uncertainties. 

Finally, we note that ref. [23] recently produced an approximated NNLO cross section by truncating 
a soft-gluon NNLL resummed calculation to order a^. Their phenomenological analysis produces 
cross sections for the LHC with extremely small scales uncertainty, of order 2-3%, sensibly smaller 
than ours. We have argued in section. [2] that such a small uncertainty also arises at the NLO+NLL 
level by requiring the factorization and renormalization scales to be equal. It will therefore be 
interesting to verify whether such reduced scale dependence found in |23j survives a test with 
independent scales, thus showing a genuine improvement due to the added NNLO terms, or whether 
it is an intrinsic consequence of keeping the scales equal. 



Note added 

After this work was completed, a new study of the ti cross section at the Tevatron and LHC 
appeared in ref. [30] . 
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A On scale and PDF uncertainties 



The definitions we have adopted in eqs. ([3]) and ^ for scale of uncertainty, and in eqs. ([7]) and ([8]) 
for PDF uncertainty, are by no means unique. In this section, we shah briefly illustrate other 
choices. 

The possibility of making different choices stems from the observation that scale and PDF uncer- 
tainties have to be combined in order to obtain an estimate of the overall uncertainty affecting the 
cross section. The way in which this combination is to be performed is at present unclear, given the 
fact that neither the scale uncertainty nor the PDF uncertainty (the latter owing to the fact that 
PDF error sets are derived in violation of the Ax^ = 1 rule) follow the laws of statistical errors. 

Scale uncertainty can in general be written as 



Act, 



Act, 



; ^(max) j.{max) 



ct(1,1), 



(30) 
(31) 



where different prescriptions can be devised for the determination of ^^(™''^^)^ ^^'""^) j q^^^^ Qf 
^^^rrtjn) ^ ^^mjn) -J ^ PDF Uncertainty is concerned, one always makes use of 

'^c''pDF+ ('^F ) Cr) - 



^ max a{setj^i) — (T{seto),a{set^i) — cj(seto), ^ 



(5(TpDF 

and then defines 



-(Cf,Cr) = (max a{seto) - a{set+i),a{seto) - a{set^i),0 ^ 



^(^PDF+ — OCrpDF+ [t,F ) SR 



AcTp 



1- / r{min) ^{min) 



(32) 
(33) 

(34) 
(35) 



where again the values of ^ ^(jnax) ^^(mm) ^ ^^™) ^ which the r.h.s. of these equations 

are evaluated are a matter of choice. 
We limit ourselves to give four examples. 

A) Our default choice, illustrated in sect. [T]and which gives rise to the results presented in this 
paper, is equivalent to solving 



max 

{€f,5r} 

min 

{€f,5r} 



a{i^,U)-a{l,l) 
c^(eF,eR)-c^(l,l) 



Uf ,Cr 



Uf S^r 



-^(1,1), 



(36) 
(37) 



{max) ^(max) 
F 1 SR 

uncertainty is defined by setting 



for (^^— j^^i'"-;) 



(mm) ^ ^(m«n) ^ ^ which are then used in eqs. ([30j) and ([3T|) . The PDF 



j-"(max) ^(max) 
SF ! ?R 



(1,1) 



j:(mm) dmin) 
?F ) ?R 



fl,l). 



(38) 
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B) The scale uncertainty is defined in the same way as done in item A ). For the PDF uncertainty, 
we set 



SF ) SR 



j-(max) j-{max) 
SF ) SR 



SF ) SR 



SF ) SR 



(39) 



with (^^^°-^) ^^^"'^)'^ a,nd ^^(^™")^ ^(™™) ^ computed again as in eqs. ([36l) and ([3 
C) We first solve 

0-(Cf,Cr) +'^<7pdf+(Cf,^r) - 



max 
{5f,Cr} 



mm 

{«f,Cr} 



/ ^(max) ^{■max)\ , c- f Amax) Amax) 

O- Uf >^R +<5'7pdf+ Uf ,^r 



0-(Cf,Cr) - '^<7pdf-(6,Cr) - 0-(l, 1) 



Uf )^R 



OCTpDF- I SF ) SR. 



^(1,1), 



^(1,1), 



(40) 



(41) 



for (^(--),^^--)) and (e^"^'"), ^i™"^ ) • These values are then used to determine the scale 
uncertainty according to eqs. ^ and and the PDF uncertainty by setting 



T{max) T{max) 
SF ) SR 



j-(max) ^{max) 
SF ) SR 



^{min) T(min) 
SF 5 SR 



j.{min) j.{min) 
Sf ) SR 



(42) 



which are then used in eqs (j34p and ([35 



D) The scale uncertainty is defined in the same way as done in item A). For PDF uncertainty, 
we solve 



max 
{?f,«r} 

max 

{«f,?r} 



I- / dmax) p{max) 

0O'poF+ I Sf ) Sr 



'^crpDF+(^F) '^R 
'^CTpDF- (^F) '^R, 

for ^ ^^ax) ^ ^^{min) ^ ^(^™*") j ^ which are then used in eqs and 



(43) 
(44) 



Items B)-D) follow the same logic, namely finding the absolute maximum and minimum of the cross 
section, by various combinations of scale and PDF uncertainties. In this sense, it is not fully justified 
to quote these two uncertainties separately, although it is still convenient for bookkeeping. These 
approaches stem from the observation that, in a hadroproduction QCD computation, unknown 
higher orders also enter the determination of the PDFs, and one is therefore entitled to use the 
full information on the PDF uncertainty in the determination of the scale dependence. In fact, the 
three methods give very similar results, with D) being the most conservative, i.e. resulting in the 
largest overall cross section uncertainty. 

On the other hand, by following the procedure outlined in item A ), one is able to better assess 
the separate dependence upon scales and PDFs. It should be observed that, while the quantities 
<5'7pDF±(?F5 ?r) depend on and roughly in the same way as the cross sections 0"(^f5?r)j the 
ratios 

5fTpDF±(CF,CR)/f^(CF,CR) (45) 
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are extremely stable with respect to variations of and ^r . This implies that the relative PDF 
uncertainty on the central value of the cross section, that is 

AapDF±/fT(l, 1) (46) 

is basically identical to any of those in eq. (j45p . if one follows item A). This is not the case for 
items B)-D); the relative uncertainty due to A(Tpdf+ (Aupdf-) tends to be larger (smaller) than 
that computed according to item A ). 

The consideration above led us to prefer the procedure of item A ) for the computation of the 
results presented in this paper. This has also the advantage that it renders the calculation less 
demanding from the point of view of CPU time. We conclude by stressing that for top production 
the procedures of items B)~D) would have given similar results as that of item A). 
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